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Abstract. The role of polar stratospheric clouds in polar ozone loss has been well 
documented. The CALIPSO satellite mission offers a new opportunity to characterize PSCs 
on spatial and temporal scales previously unavailable. A PSC detection algorithm based on a 
single wavelength threshold approach has been developed for CALIPSO. The method 
appears to accurately detect PSCs of all opacities, including tenuous clouds, with a very low 
rate of false positives and few missed clouds. We applied the algorithm to CALIPSO data 
acquired during the 2006 Antarctic winter season from 13 June through 31 October. The 
spatial and temporal distribution of CALIPSO PSC observations is illustrated with weekly 
maps of PSC occurrence. The evolution of the 2006 PSC season is depicted by time series of 
daily PSC frequency as a function of altitude. Comparisons with ‘virtual’ solar occultation 
data indicate that CALIPSO provides a different view of the PSC season than attained with 
previous solar occultation satellites. Measurement-based time series of PSC areal coverage 
and vertically-integrated PSC volume are computed from the CALIPSO data. The observed 
area covered with PSCs is significantly smaller than would be inferred from a temperature- 
based proxy such as T NA t but is similar in magnitude to that inferred from Tsts. The potential 
of CALIPSO measurements for investigating PSC microphysics is illustrated using 
combinations of lidar backscatter coefficient and volume depolarization to infer composition 
for two CALIPSO PSC scenes. 


1 Introduction 

Extensive measurements and modeling activity over the past two decades have firmly 
established that polar stratospheric clouds (PSCs) play an essential part in the springtime 
chemical depletion of ozone at high latitudes (Solomon, 1999). The role of PSCs is twofold. 
First of all, PSC particles provide the catalytic surfaces for heterogeneous chemical reactions 
that transform stable chlorine and bromine reservoir species into highly reactive ozone- 
destructive forms. Secondly, if PSC particles grow sufficiently large to sediment, they can 
irreversibly redistribute gaseous odd nitrogen (a process commonly known as denitrification), 
thereby slowing the reformation of the chlorine reservoirs and prolonging the ozone depletion 
process. In spite of the recent progress made in understanding PSCs and their chemical 
effects, there are still significant gaps in our knowledge, primarily in the areas of solid particle 
formation and their denitrification potential. These uncertainties limit our ability to 
accurately represent PSCs in global models and call into question the accuracy of PSC 
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predictions for a future stratosphere whose composition differs from the current state. This is 
of particular concern in the Arctic, where winter temperatures hover near the PSC threshold 
and, hence, future stratospheric cooling could lead to enhanced cloud formation and 
substantially greater ozone losses. 


Much of our present understanding of PSCs is based on remote sensing observations. Long- 
term solar occultation data records from the Stratospheric Aerosol Measurement (SAM) II 
(Poole and Pitts, 1994) and Polar Ozone and Aerosol Measurement (POAM) II and III 
(Fromm et al., 1997; Fromm et ah, 1999; Fromm et ah, 2003) instruments have established 
the general climatology of PSCs such as their occurrence at very cold temperatures (T < 195- 
200 K) and their spatial distribution and seasonal variability in both the Arctic and Antarctic. 
Measurements from polarization sensitive airborne and ground-based lidar systems have also 
made valuable contributions to our understanding of PSCs. In particular, lidar observations 
have established that there are three primary forms of PSC particles: supercooled liquid 
ternary (H 2 SO 4 -HNO 3 -H 2 O) solution droplets (STS) that evolve from background aerosol 
without a phase change; solid nitric acid hydrate crystals, most likely nitric acid trihydrate 
(NAT); and FLO ice (Poole and McCormick, 1988; Browell et ah, 1990; Toon et ah, 1990). 
Multi-year ground-based lidar observations have been utilized to study the seasonal 
chronology of PSC composition (David et al, 1998; Gobbi et ah, 1998; Santacesaria et al. 
2001; Adriani et ah, 2004). More recently, satellite limb emission measurements of the 
infrared spectra of PSCs by the Michelson Interferometer for Passive Atmospheric Sounding 
(MIPAS) instrument on ENVISAT have provided valuable new insight to PSC composition 
(e.g., Spang et ah, 2005; Hopfner et al, 2006). Both solar occultation and lidar data, however, 
have inherent shortcomings that limit a more comprehensive understanding of PSC processes. 
For instance, solar occultation measurements have very coarse spatial resolution (hundred of 
kilometers), are limited to 14-15 profiles per day in each hemisphere, and only occur at the 
day-night terminator. Ground-based lidar records are of course from single locations and are 
often interrupted by optically-thick tropospheric clouds, while airborne lidar PSC data records 
are both short-term and limited to relatively short along-track spatial coverage. 


Spacebome lidar measurements are particularly well-suited for probing PSCs and may 
provide a more comprehensive picture of clouds than attainable with these existing PSC 
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remote sensing data sources. The Lidar-In-space Technology Experiment (LITE) 
demonstrated the potential of space lidar for observation of clouds and aerosol (Winker et al., 
1996). The Geoscience Laser Altimeter System (GLAS), launched aboard the Ice Cloud and 
land Elevation Satellite (ICESat) in 2003, measures atmospheric backscatter profiles at two 
wavelengths (532 and 1064 nm) and has observed PSCs (Palm et ah, 2005), but this system is 
not polarization sensitive and due to laser lifetime concerns the observational periods are 
limited to four 33-day periods per year (Abshire et ah, 2005). The Cloud-Aerosol-Lidar and 
Infrared Pathfinder Satellite Observations (CALIPSO) mission, launched on 28 April 2006, 
provides a unique set of measurements to study PSCs. The CALIPSO PSC data record will 
have much higher vertical and spatial resolution than solar occultation data, but more 
importantly, the CALIPSO data are collected continuously along the orbit track so that PSC 
information will be extracted in the polar night region and not solely at the terminator. 
Similar to many airborne and ground-based systems, the CALIPSO backscatter data are 
recorded in two orthogonal polarization channels at 532 nm, and in a third channel at 
1064 nm. This combination of measurements provided by CALIPSO offers promise of 
additional information on microphysical properties, particularly the presence of large solid 
particles in mixed-phase clouds, a crucial element for understanding and accurate prediction 
of the role that PSCs play in polar ozone loss, but without interference from tropospheric 
clouds and the limits in spatial and temporal coverage of ground-based and airborne lidar data 
records. 


In this paper, we describe an algorithm to detect polar stratospheric clouds in the CALIPSO 
data. We then apply the algorithm to CALIPSO data from the 2006 Antarctic winter season 
to demonstrate the robustness of the approach and highlight the unique capabilities of 
CALIPSO for characterizing PSCs. 


2 CALIPSO Description 

The CALIPSO mission is designed to provide global, vertically resolved measurements of 
clouds and aerosols to improve our understanding of their role in climate forcing (Winker et 
al., 2003). The CALIPSO instrument payload consists of the Cloud-Aerosol Lidar with 
Orthogonal Polarization (CALIOP), the Infrared Imaging Radiometer (HR), and the Wide 
Lield Camera (WLC). CALIOP is a two-wavelength, polarization sensitive lidar that provides 
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high vertical resolution profiles of backscatter coefficient at 532 and 1064 nm, as well as two 
orthogonal (parallel and perpendicular) polarization components at 532 nm (Winker et al., 
2007). Table 1 provides a summary of the CALIOP instrument parameters. Although the 
fundamental sampling resolution of CALIOP is 30 m in the vertical and 333 m in the 
horizontal, an altitude-dependent on-board averaging scheme is employed that provides 
highest resolution in the troposphere and lower resolution in the stratosphere. The resultant 
downlink resolution of the CALIOP data is listed in Table 2. A detailed discussion of the 
CALIOP data products can be found in Vaughan et al. (2004). 


CALIPSO is in a 98° inclination orbit with an altitude of 705 km and is a part of the Aqua 
satellite constellation that includes the Aqua, CloudSat, Aura, and PARASOL platforms. 
Although PSCs are not one of its primary targets, CALIPSO is an ideal platform for studying 
polar processes since its orbit provides extensive measurement coverage over the polar 
regions of both hemispheres with an average of fourteen orbits per day extending to near 82° 
latitude. Figure 1 illustrates the measurement coverage over the Antarctic for a single day (28 
June 2006) with blue points representing the nighttime measurement locations and red points 
representing the daytime measurement locations. On average, over 300,000 lidar profdes are 
acquired per day at latitudes poleward of 55°, providing a unique dataset for studying the 
occurrence, composition, and evolution of PSCs. In contrast, SAM II made less than 60,000 
observations over the Antarctic during the entire eleven year period from 1978 to 1989. 


3 PSC Detection Algorithm 

The PSC detection algorithm is an automated process to identify PSCs in the CALIOP Level 
IB Data Product. In the initial version described herein, we have adopted a single wavelength 
threshold approach using the CALIOP 532-nm backscatter coefficient measurements although 
we anticipate future versions will also incorporate the CALIOP 532-nm depolarization and 
1 064-nm backscatter coefficient measurements. All results shown in this paper are based on 
Version 1.10 of the CALIPSO data product available publicly through the NASA Langley 
Atmospheric Science Data Center (http://eosweb.larc.nasa.gov/) . 
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For this initial PSC detection algorithm, we restricted our analysis to CALIOP nighttime 
measurements only. While PSCs can be seen in the daytime data, we find that the daytime 
measurements have higher background signals due to scattered sunlight. This forces the use of 
a significantly higher PSC detection threshold for the CALIOP daytime measurements than 
for the nighttime measurements. We will consider the use of daytime measurements in the 
future. Figure 2 shows the distribution of the nighttime measurement locations as a function 
of latitude for each month of the PSC season. As expected for CALIPSO’s orbit inclination 
of 98°, during most of the season there is a sharp maximum in the number of measurements at 
the crossover point of the orbit near 82° S with the remainder of the measurements more 
evenly distributed in latitude. By October, the Sun is illuminating a larger portion of the polar 
region and the majority of nighttime measurement locations are distributed between about 55° 
and 75° S. The latitudinal distribution of nighttime measurement locations will be taken into 
consideration when we examine the statistics of PSC occurrence over the Antarctic region as a 
whole. 


3.1 Data Preparation 

The standard CALIOP Level IB data files contain a half orbit (day or night) of calibrated and 
geolocated full resolution profiles of attenuated backscatter coefficient at 532 and 1064 nm 
and perpendicular attenuated backscatter coefficient at 532 nm. In addition, the Level IB files 
contain profiles of meteorological data that are also components of the detection algorithm. 
This data is derived from the Goddard EOS Data Assimilation System (GEOS-4) products 
(Bloom et al., 2005; Lin, 2004) interpolated to the CALIPSO measurement locations. 


The altitude region of interest for PSC detection ranges from near the tropopause up to about 
30 km and spans several CALIOP spatial resolution regimes (see Table 2). Specifically, the 
resolution of the Level IB CALIOP data is 60-m vertical and 1-km horizontal for altitudes 
between 8.2 km and 20.2 km and 180-m vertical and 1 .674cm horizontal for altitudes between 
20.2 km and 30 km. This change in averaging scales produces a distinct difference in the 
noise characteristics of the data across the 20.24cm boundary as illustrated in Figure 3a. In 
addition to having finer spatial resolution, the data below 20.2 km also exhibit a significantly 
larger range in background noise. Since the cloud detection algorithm is based on ensemble 
statistics, mixing data with significantly different noise characteristics would produce 
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unpredictable results. One solution would be to treat the two altitude regimes separately for 
cloud detection, but this would add significant complexity to the algorithm and it is possible 
that the two methods would give different results and create an artificial boundary in the PSC 
analyses. Instead, we applied additional spatial averaging to all data below 20.2 km to closely 
match the resolution of the data above 20.2 km. The additional averaging produces a dataset 
with consistent resolution and noise characteristics over the entire altitude range of interest as 
shown in Figure 3b. The change in spatial scales and the concomitant change in noise 
character at the 8.2, 20.2, 30.1 km boundaries are sometimes visible in the CALIPSO lidar 
Level IB browse images (available at http://www-calipso.larc.nasa.gov/products/lidar/) . 


Since the lower bound for PSC opacity merges into the background aerosol, we have chosen 
to perform additional averaging of the data to reduce the noise levels in the data and lower the 
detection limit for PSCs. Some of this noise occurs in the form of radiation induced noise 
spikes. Although these events are typically isolated, they occur with sufficient frequency that 
their inclusion in our analyses would produce systematic bias in the ensemble statistics and 
false positive cloud detection. We found that first applying a 5-km horizontal median filter 
eliminates the majority of these noise spikes while only modestly impacting the spatial 
resolution. To further reduce the noise, we also averaged the data to 540-m resolution in the 
vertical. 


We applied the additional processing described above to the Level IB CALIOP data to 
produce an input dataset for the cloud detection algorithm that has single profile resolution of 
540-m vertical and 5-km horizontal. To reduce the data volume, we limited our analyses to 
latitudes poleward of 55° S and to altitudes between 10 and 30 km. 


3.2 Algorithm Description 

The CALIPSO PSC detection algorithm is similar to the traditional single wavelength 
extinction threshold techniques used in solar occultation PSC identification schemes such as 
Poole and Pitts (1994) and Fromm et al. (2003). In these previous studies, PSCs were 
identified as those measurements having extinction coefficients significantly larger than the 
background (non-PSC) aerosol extinction. Instead of extinction, the CALIPSO cloud 
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detection algorithm is based on measurements of scattering ratio at 532 nm, R(532), defined 


as 


£(532) - 
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where /3 m (532) is the volume molecular backscatter coefficient at 532 nm and /3t(532) is the 
total (molecular + aerosol) volume backscatter coefficient at 532 nm. J3 m (532) is calculated 
from the GEOS-4 molecular number density profiles provided in the CALIOP Level IB data 
files and a theoretical value for the molecular scattering cross section as described in Hostetler 
et al. (2006). (3t(532 ) is derived directly from the attenuated backscatter coefficient at 532 nm 
reported in the CALIOP Level 1 B data files by applying a first-order correction for molecular 
and ozone attenuation based on GEOS-4 molecular and ozone number density profiles 
provided in the CALIOP Level IB data files. 


Figure 4 shows an ensemble of CALIOP scattering ratio measurements as a function of 
observed temperature from 13 June 2006 over the Southern polar region. Within the winter 
polar vortex, the background aerosol population is characterized by very low optical depth 
values (Thomason and Poole, 1993) and is represented in Figure 4 by the nearly vertical 
family of points with scattering ratio values close to 1 . The spread of the background aerosol 
data is representative of the magnitude of the noise in the measurements (g= 0.32) at these 
spatial scales (540-m vertical x 5-km horizontal) and is consistent with the observed CALIOP 
nighttime signal-to-noise ratio in the stratosphere. The statistics of the data totally change 
below about 195 K where the upper limit of the scattering ratio values increases dramatically 
from ~2 to ~8. In Figure 4, the enhanced scattering ratio values in excess of ~2 observed at 
cold temperatures roughly correspond to the PSC visible in Figure 3. This characteristic 
dependence of PSC occurrence on cold temperature is the physical basis for our detection 
algorithm. It allows us to characterize the statistics of the background aerosol data at warmer 
temperatures apart from the presence of PSCs. PSCs are then identified as statistical outliers 
in terms of scattering ratio from the background aerosol ensemble occurring at temperatures 
below a pre-selected threshold temperature, T PS c. The background aerosol population is 
characterized by the ensemble of points that occur at temperatures warmer than Tpsc, with the 
assumption that the optical properties of the background aerosol are not significantly 
dependent on temperature. Based on the distribution of the background aerosol population, a 
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scattering ratio threshold, Rpsc, is defined such that a large fraction of the warm temperature 
points lie below this value. Points that exceed this value and are below Tpsc are considered to 
be PSCs. The ensemble statistics and corresponding PSC scattering ratio threshold, Rp, are 
computed on a complete day of nighttime CALIOP observations, then PSCs are identified on 
a profile-by-profile basis. 


The algorithm is summarized by the following steps: 

1 . The background aerosol population is defined as the ensemble of scattering 
ratio points with temperatures above Tpsc =198 K. 

2. Rx is defined as the scattering ratio value corresponding to the 99.5 percentile 
(~ 2.6 a) of the background aerosol population defined in step (1). 

3. PSCs as identified as those points with temperature below Tpsc and scattering 
ratio values exceeding R T . 


The detection process is illustrated graphically in Figure 5a for the same ensemble of 
scattering ratio points as shown in Figure 4. The scattering ratio distribution for the entire 
ensemble is represented by the solid black line and is characterized by a roughly Gaussian- 
shaped core corresponding to the background aerosol and a long positive tail corresponding to 
the PSCs. The background aerosol distribution alone, defined by all points with T > Tpsc, is 
depicted by the black dotted line and is nearly Gaussian shaped with no positive tail and a 
mode scattering ratio nearly identical to that of the entire ensemble. The dashed red line 
indicates the value of Rt corresponding to the 99.5 percentile of the background aerosol 
distribution alone. The 99.5 percentile value was chosen as the PSC threshold based on trial 
and error. The sensitivity of the detection algorithm to both Rx and Tpsc are discussed in 
Section 3.3. The blue line represents the distribution of ensemble points with T < T PS c, which 
includes both background aerosol and PSCs. There is a peak in the cold (blue) distribution 
corresponding to the background aerosol at cold temperatures and a pronounced positive tail 
corresponding to the PSCs. The slight positive shift of the mode scattering ratio of the cold 
background points (blue line) is suggestive of growth of the entire aerosol population as the 
temperature cools. This may be a slight violation of our basic assumption that the ensemble 
optical properties of the background aerosol are not significantly dependent on temperature. 
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However the shift is small and we are forced to live with this behavior. By definition, all 
points in the cold (blue) distribution with scattering ratios larger than Rt are identified as 
PSCs. For contrast, the corresponding ensemble distribution curves for a day in late October 
with no PSCs present are shown in Figure 5b. In the absence of PSCs, all three curves are 
representative of the background aerosol distribution and have similar near-Gaussian shapes 
with no positive tails. 


As is evident from Figure 4, PSCs with scattering ratios greater than about 2 are easy to 
discriminate from the background aerosol population. However, optically thin PSCs with 
backscatter ratios less than about 2 are more difficult to detect above the measurement noise 
present in the 5 -km resolution backscatter data. To enhance the detection of these more 
tenuous PSCs, two additional passes are made through the data with 25-km and 75-km 
horizontal averaging applied to the data. Each successive level of increased spatial averaging 
reduces the noise level in the data, exposing the more tenuous PSCs. Figure 6 illustrates the 
effect of the successive averaging on the ensemble distributions and scatter plots. Compared 
with Figures 4 and 5 a, each level of additional smoothing reduces the spread in the data, 
lowering the scattering ratio threshold for PSC detection. On average, Rt has a value of near 
1.8 for the 5-km resolution data; 25-km averaging reduces R T to about 1.3; and 75-km 
smoothing further reduces Rj to values near 1.2. Using this successive averaging approach, 
the threshold for PSC detection is decreased by almost a factor of four. 


3.3 PSC Detection Sensitivity 

The detection algorithm is empirical by nature and the selection of values for Rj and Tpsc is 
not clear cut. Since the value of Tp§c defines the ensemble of scattering ratio points used to 
characterize the optical properties of the background aerosol population, the most important 
criterion for selecting Tpsc is to choose a value sufficiently large to exclude PSCs from the 
background aerosol ensemble. We chose a value of T PS c = 198 K. This value is warm enough 
to exclude most if not all clouds from the background ensemble. We found that changing 
Tpsc by +1-2° changes the total number of PSCs identified by about 1%. However, dropping 
the PSC temperature threshold to 196 K clearly allows some PSC values to be included in the 
background aerosol ensemble. As a result, this raises the PSC scattering ratio threshold and 
decreases the number of PSCs identified by as much as 6%. It is crucial to set the threshold 
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temperature large enough to exclude most if not all PSCs while not setting it so large that the 
sample size becomes too small. We have found that 198 K is sufficient in all cases. The 
identification of PSCs is also sensitive to the percentile value selected for the scattering ratio 
threshold. We considered using the 97.5 percentile (~2 a) and the 99.85 percentile (~3 a), as 
well as, the 99.5 percentile value as the scattering ratio threshold. Visual inspection of the 
data revealed that using the 97.5 percentile produced a large number of false positives. These 
appeared in the images as isolated PSC ‘pixels’ highly indicative of noise rather than the 
coherent structure of a geophysical object. On the other hand, the 99.85 percentile results 
were similar to the 99.5 percentile, but reduced the overall size of the sample to the point that 
determining the threshold Rt value became noisy and highly variable with time. Therefore, 
we find that the 99.5 percentile definition is robust in the sense that it produces few false 
positives and allows a sample size sufficiently large to produce robust statistics. 


Two examples are shown to illustrate the performance of the PSC detection algorithm. 
Figure 7a shows the CALIOP 532-nm backscatter coefficient data from a single orbit on 15 
June 2006. The corresponding PSC mask produced from the PSC detection algorithm is 
shown in Figure 7b. The averaging scale (5-, 25-, or 75-km) at which the PSC was detected is 
indicated by the shading. For this optically-thin PSC case, the detection algorithm is able to 
discriminate even the thinnest parts of the cloud from the background signal, although higher 
levels of spatial averaging are required to detect the more tenuous regions above the noise. 
Figure 8 show a similar comparison for an optically thick PSC case from 24 July 2006. As 
expected, the optically-thick PSCs are easily identified by the detection algorithm during the 
first pass at the 5 -km resolution. It is noteworthy that the cloud edges are accurately 
reproduced and the tenuous clouds surrounding the dense cloud are also being detected during 
subsequent passes at larger spatial scales. Based on visual inspection of dozens of images 
similar to those shown in Figures 7 and 8, the CALIPSO PSC detection algorithm appears to 
be working exceptionally well. Based on the frequency of random, isolated PSC points, the 
rate of false positives is extremely low, 0.1% or less. Missed clouds are more difficult to 
quantify, but following visual inspection of a number of images we did not observe any 
features likely to be PSCs that weren’t identified as such. At the lower end of the scale, PSC 
scattering ratio merges into that of the background aerosol. While higher levels of smoothing 
capture much of this behavior, we are, by definition, missing some segment of this low 
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backscatter continuum. In addition, PSCs that are characterized by very low particle number 
densities such as a NAT ‘haze’ may not produce detectable enhancements in backscatter 
coefficient and could therefore also be missed in this initial version of the algorithm. Future 
versions will attempt to reduce this aerosol-PSC ambiguity by incorporating the additional 
1064-nm backscatter coefficient and polarization information available in the CALIPSO data 
products. For instance, inclusion of the 532-nm perpendicular backscatter coefficient data 
may allow us to detect thin layers of NAT haze. 


4 2006 Antarctic PSC Season 

As a test of the PSC detection algorithm, we examined the CALIPSO PSC observations 
during the 2006 Antarctic season. We applied the PSC detection algorithm to all CALIPSO 
data from 13 June (first date available) until 31 October 2006. To quantify the occurrence of 
PSCs spatially, we defined a 21 x 21 -box grid in polar stereographic coordinates centered on 
the South Pole and extending to 55° S at its sides. Each of the 441 grid boxes covers an area 
approximately 350 km x 350 km. For a selected time interval such as a day or week, a PSC 
frequency is defined for each grid box as the ratio of the number of PSCs observations to the 
total number of CALIPSO observations within that box. Later, we will utilize this PSC 
frequency as roughly the equivalent to areal coverage of PSCs within a grid box. The 
resolution of the grid was chosen as a compromise between the need for adequate sampling 
within each grid box and the desire to retain high spatial resolution. This is more of an issue 
for the low latitude boxes. For example, in a typical week, the grid boxes at high latitudes 
have nearly two thousand CALIOP observations, while grid boxes at lower latitudes near the 
edge of the grid have several hundred observations. 


4.1 Seasonal Evolution 

The meteorological conditions in the Antarctic stratospheric during the 2006 winter set the 
stage for a record ozone hole in terms of both area and ozone mass deficit (Gutro, 2006; 
WMO Bulletin, 2007). Temperatures poleward of 50° S at 50 mb were colder than the 1979- 
2005 average from early August through November. During this same time period, the area 
of the polar vortex at 450 K was significantly larger than average. This led to a much larger 
than average ozone hole and, if fact, during many days in November the ozone hole area was 
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larger than had been observed at any time between 1979 and 2005 (WMO Bulletin, 2007). In 
addition, the ozone mass deficit was higher than for any year since 1997 (WMO Bulletin, 
2007). Since CALIPSO began routine data acquisition on 13 June, the initial onset of PSCs 
was not observed by CALIPSO in 2006. 


To examine the spatial distribution of PSCs over the course of the 2006 Antarctic winter, 
weekly PSC frequency maps were produced for every 540 m in altitude between 10 and 
30 km. Examples from two altitudes are shown in Figures 9 and 10. At 20 km (Figure 9), 
PSCs were prevalent over a large fraction of the polar region from the onset of the CAFIPSO 
measurements in mid- June until mid-August, after which PSCs significantly diminished and 
eventually disappeared by mid-October. At 15 km (Figure 10), PSCs were less widespread 
early in the season than at 20 km, but became ubiquitous over the polar region by July and 
persisted well into October. At all altitudes, PSCs were commonly observed over most of the 
Antarctic continent, but their frequency rapidly diminished away from the continent at lower 
latitudes. Although the location of the maximum PSC frequency tended to vary from week to 
week, there were clearly some preferred regions for PSC occurrence. This is evident in 
Figure 11, where we show the longitudinal distribution of PSCs as a function of altitude. 
PSCs were observed over all longitudes, but the two favored regions for PSC occurrence were 
over the ice sheets of East Antarctica and over the Antarctic Peninsula near 90° W where the 
PSC frequency occasionally exceeded 60%. The longitudinal distribution of temperature 
(from GEOS-4) confirms that the longitudinal distribution of PSCs is strongly correlated with 
cold temperatures as shown by Poole and Pitts (1994). 


To examine the temporal evolution of PSCs for the Antarctic as a whole, we produced daily 
time series of PSC frequency as a function of altitude. In this context, the PSC frequency 
represents the cumulative sum over the entire grid of the PSC frequency (areal coverage) 
within each grid box on a given day at a given altitude. To account for both the uneven 
latitudinal sampling (see Figure 2) and the tendency for PSCs to occur more frequently at 
higher latitudes, the PSC frequencies are weighted by their location in either the inner core or 
outer fringe of the grid. The inner core consists of the central 11x11 grid boxes and the outer 
fringe contains all remaining grid boxes. The calculated PSC statistics for these two regions 
are given weights proportional to their relative areal size: 121/441 for the inner core and 
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320/441 for the outer fringe. The CALIPSO time series of daily PSC frequency is shown in 
Figure 12a and the corresponding daily-averaged GEOS -4 temperature at the CALIPSO 
measurement locations is shown in Figure 12b. A maximum PSC frequency of slightly less 
than 0.30 (interpreted as 30% of the grid was covered with PSCs) occurred during several 
periods between the end of June and mid-August. By early September, PSC frequency had 
diminished significantly at all altitudes. Then in mid-September, PSC activity increased again 
for a few weeks before tapering off in October. The altitude of the peak PSC occurrence 
slowly decreased from near 22 km in June to below 15 km by September and October. As 
expected, the maximum PSC frequency is highly correlated with region of coldest 
temperatures. 


4.2 Comparison with Solar Occultation PSC Climatologies 

Much of our current knowledge of PSC occurrence over the Antarctic has been based on the 
historical records of solar occultation satellite instruments such as SAM II (1978-1993) and 
POAM II/II data (1996-2006) (e.g., Poole and Pitts, 1994; Fromm et al., 1997; Nedoluha et 
al., 2003). A significant limitation of solar occultation climatologies is that they are 
representative of only the temporally varying latitudes at which the occultation observations 
are being made and not of the entire polar region, especially the polar night which is never 
sampled by the solar instruments. In general, solar occultation measurement latitudes for sun- 
synchronous platforms like POAM III and SAM II vary slowly (nominally l°-2° per week) 
with season, reaching highest (lowest) latitudes near the equinoxes (solstices). Solar 
occultation instruments will generally saturate when observing optically thick PSCs and, 
hence, this class of clouds is not included in their data records. In contrast, no saturation is 
apparent in any of the CALIPSO PSC images. With these limitations, there have been 
questions regarding how generally representative these solar occultation climatologies are. 
With CALIPSO’s high-density spatial sampling, it is straightforward to mimic solar 
occultation sampling and produce ‘solar occultation’ -like time series of PSC frequency. 
Using SAM II measurement latitudes as the reference for solar occultation, we produced a 
time series of daily PSC frequency using only CALIPSO measurements from the latitude 
band that would have been observed by SAM II. This virtual ‘solar occultation’ PSC 
frequency time series is shown in Figure 13a. Compared to the Poole and Pitts (1994) 
climatology for the Antarctic, the virtual ‘solar occultation’ time series is similar with peak 
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frequency occurring below 16 km in September. In fact, the 2006 season is very similar to the 
cold 1987 season statistics shown in Poole and Pitts (1994). This is quite different than the 
CALIPSO time series (Figure 12) where with the peak PSC frequency occurs much earlier 
near the beginning of August and higher in altitude near 20 km. We also constructed a virtual 
‘solar occupation’ time series representing the POAM II/III measurement latitudes, but both 
instruments sampled at higher latitudes in September than observed by CALIPSO so a 
complete POAM-like time series was not possible. The virtual POAM solar occultation time 
series (excluding September) is similar to the evolution observed by POAM II/III as shown in 
Fromm et al. (1997) and Nedoluha et al. (2003) with peak PSC frequency observed later in the 
season than observed with CALIPSO. Considering the solar occultation sampling pattern, it 
is no surprise that their observed maximum in PSC frequency occurs in September since that 
is when the occultation measurement locations are near their highest latitude and the observed 
temperatures are coldest (Figure 13b). In light of the differences between the CALIPSO and 
virtual ‘solar occultation’ PSC statistics, it is clear that in 2006, at least, the temporal 
distribution of solar occultation PSC observations would not be representative of the Southern 
polar region as a whole. 


4.3 PSC Processing Diagnostics 

Chemical ozone loss in the polar regions is strongly related to the geographical extent of PSCs 
during the winter (e.g., Rex et al., 1999, 2004, 2006; Tilmes et al., 2004). To date, realistic 
measurement-based estimates of PSC areal coverage have not been feasible due in large part 
to the limited sampling pattern of solar occultation instruments. Instead, proxies for total PSC 
area have been developed based on temperature. The most frequently used proxy for PSC 
areal extent is the area with temperature less than Tnat calculated from global temperature 
analyses such as ECMWF or GEOS-4 and an assumption of thermodynamic equilibrium with 
NAT, the most stable phase of PSCs that occur at temperatures above the water frost point 
(Hanson and Mauersberger, 1988). However, with the extensive measurement coverage 
afforded by CALIPSO, observational-based estimates of the geographical extent of PSCs are 
now possible. The daily PSC frequency (Figure 12a) easily scales into total PSC areal 
coverage as a function of altitude and is shown in Figure 14. By early July, PSCs are 
covering over 10 million km 2 at altitudes between about 20-23 km. The region of maximum 
PSC area expands in size into early August when PSCs are covering over 10 million km 2 at 
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altitudes between about 15 and 23 km. The total volume of atmosphere encompassed by 
PSCs can be determined by vertical integration of the PSC area and is shown in the bottom 
panel of Figure 14. The PSC volume generally increases from mid-June until it reaches a 
peak of near 150 million km 3 on about 1 August, after which it steadily decreases until 1 
September. There is a slight increase in PSC volume during the last two weeks of September, 
corresponding to the brief upturn in PSC frequency at that time, after which PSC volume 
continues to decrease until it reaches zero by the end of October. 


The CALIPSO-based estimates of PSC areal coverage provide a means to evaluate the 
robustness of temperature proxies for PSC formation such as Tnat- We used the formulation 
of Hanson and Mauersberger (1988) along with profiles of HNO 3 and H 2 O mixing ratio from 
the EOS (Aura) Microwave Limb Sounder (MLS) (Data Version 1.5) (Froidevaux et al., 
2006) to produce maps of Tnat- To account for the large spatial and temporal variations in 
the observed HNO 3 and H 2 O mixing ratios, we produced daily-averaged, zonal mean profiles 
of HNO 3 and H 2 O mixing ratio for 5° latitude bands from 55° S to 85° S. Figures 15 and 16 
show sample time series of Aura MLS zonal HNO 3 mixing ratio, H 2 O mixing ratio, and the 
derived Tnat for the 5° latitude bands centered at 67.5° S and 82.5° S. During the 2006 
winter, Aura MLS observed very low HNO 3 mixing ratios (< 2 ppbv) and H 2 O mixing ratios 
(< 3 ppmv) over a large area within the polar vortex. At 82.5° S, significant denitrification 
and dehydration had occurred by mid-winter and, as a result, T N at was below 190 K over a 
large altitude region during much of August and September. At 67.5° S, the onset of 
denitrification was delayed and not nearly as severe as what occurred at higher latitudes and 
the corresponding Tnat was greater than 190 K over most of the stratosphere. Using the MLS- 
derived zonal mean profiles of Tnat, an estimate of the area with T < Tnat was produced in 
an analogous fashion as the CALIPSO PSC area by calculating the cumulative sum over the 
entire grid of the fraction of occurrence within each grid box of GEOS-4 temperatures less 
than Tnat- The result is shown in Figure 17. Qualitatively, the temporal and spatial evolution 
of the area with T < Tnat is similar to the observed CALIPSO PSC area, but its magnitude is 
larger than the CALIPSO PSC area by almost a factor of two. The accuracy of the Tnat area 
calculations are dependent on the accuracy of the GEOS-4 temperature and MLS gas species 
mixing ratio data, but it is unlikely these uncertainties are large enough to explain the 
observed discrepancy. It is also possible that the CALIPSO PSC areal extent is being 
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underestimated since a subset of optically-thin PSCs such as NAT haze may not be detected 
in this version of the algorithm. 


An alternative proxy for PSC occurrence is the equilibrium temperature of STS, Tsts- We 
produced an estimate of the area with T < Tsts (shown in Figure 18) where we have assumed 
that Tsts is 4 K colder than Tnat (Tabazadeh et al., 1994). While there are still differences, 
the evolution of the area with T < Tsts agrees much better with the CALIPSO PSC area and is 
nearly identical in magnitude. Therefore, for the 2006 Antarctic winter, there is much better 
agreement between the observed CALIPSO PSC area and Tsts than Tnat- This is in 
approximate agreement with Santee et al. (1998) who showed strong correspondence between 
the area of T < 192 K and the area of gas-phase HNO3 loss, an indicator of PSC formation, 
but weak correspondence between the area of T < 195 K and the area of gas-phase HNO3 loss. 
The CALIPSO-based PSC areas will provide an opportunity to further investigate the 
accuracy of these proxies and possibly develop improved predictors of PSC occurrence. 


4.4 PSC Composition Studies 

Although not a primary goal of this work, the CALIPSO lidar measurements are also well- 
suited for probing the composition of PSCs. The CALIOP backscatter data are recorded in 
two orthogonal polarization channels at 532 nm, and in a third channel at 1064 nm. This 
combination of data provides crucial information of PSC particle size (through changes in 
absolute backscatter and its wavelength dependence) and phase (via the 532-nm perpendicular 
backscatter. For instance, the presence of large non-spherical particles in sufficient quantities 
will produce changes in the ensemble optical properties such as the ratio of 1064-nm to the 
532-nm backscatter (color ratio) and especially the 532-nm depolarization ratio. For 
illustrative purposes, we show a few examples of how the CALIPSO data may be exploited to 
infer PSC composition. In a follow-up paper by (Poole et al., 2007), we explore the PSC 
microphysical information content of the CALIOP data in detail. 


In Figure 19 is shown a scatter diagram of CALIOP scattering ratio at 532 nm versus volume 
depolarization at 532 nm for measurements identified as PSCs on a single orbit from 24 July 
2006. A commonly used scheme for optical classification of PSC remote sensing 
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observations has evolved from analyses of polarization sensitive lidar data (Poole and 
McCormick, 1988; Browell et al, 1990; Toon et al., 1990): type la for NAT, type lb for STS, 
and type 2 for H 2 O ice. More recent studies by Toon et al. (2000) and Biele et al. (2001) have 
shown that such “typing” of PSC data is somewhat of an oversimplification, with the clouds 
often being mixtures of solid and liquid particles in reality. For illustrative purpose, we subset 
the CALIPSO PSC data into three groups (STS, FFO ice, and mixtures of STS, FLO ice, and 
NAT) based on the combinations of backscatter ratio and volume depolarization as shown in 
Table 3. The CALIPSO measurements roughly corresponding to STS (red) and FLO ice 
(blue) are indicated in the figure. All other points are considered mixtures. Figure 20 shows 
two examples of PSC composition inferred from the CALIPSO measurements. The image in 
Figure 20a corresponds to the tenuous PSC shown in Figure 7. Based on this simple 
classification scheme, the PSC consists almost exclusively of STS droplets. The image in 
Figure 20b corresponds to the intense PSC event shown in Figure 8. In this case, the cloud 
consists of a large FLO ice core that is surrounded by mixtures of particles. Near the bottom 
right edge of the cloud, a region of mostly STS is observed. This particular PSC occurred near 
the Antarctic Peninsula and is likely a mountain-wave induced cloud. Although we have 
used a simplistic analysis in these examples, it is clear that the CALIPSO data offers promise 
of additional information on PSC microphysical properties. 


5 Summary and Conclusions 

The CALIPSO mission provides a unique dataset for studying the formation, evolution, and 
composition of PSCs. The CALIPSO PSC data record has higher vertical and spatial 
resolution than previously attainable from spaceborne remote sensors and extends into the 
polar night region. As the first step in exploiting this promising database, we have developed 
an automated algorithm to detect PSCs in the CALIPSO Level IB data. The detection 
approach is based on CALIOP scattering ratio measurements at 532 nm and uses a single 
threshold in scattering ratio for PSC identification. Successive levels of smoothing are 
applied to the data to increase the sensitivity for detecting optically thin PSCs. Based on 
visual inspection of dozens of CALIPSO PSC cloud images, the detection algorithm is robust 
with only a small fraction of false positives and few if any missed clouds that were apparent 
in the backscatter data. 
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As a test of the algorithm, we examined CALIPSO PSC observations from the 2006 Antarctic 
winter (13 June until 31 October). With its extensive measurement coverage in the polar 
regions, CALIPSO provides a comprehensive view of the spatial and temporal distribution of 
PSCs during the Antarctic winter. CALIPSO observed widespread PSCs over much of the 
Antarctic region from June through August, after which PSCs significantly diminished, except 
for a brief resurgence in the latter half of September, and eventually disappeared by mid- 
October. Throughout the season, PSCs were most frequently observed over East Antarctica 
and near the Antarctic Peninsula, where at times the PSC frequency exceeded 60%. 
Comparison of the CALIPSO PSC statistics with virtual ‘solar occupation’ data indicates that 
solar occultation PSC statistics would not be representative of the Antarctic region as a whole 
during the 2006 winter. Although this is not surprising given their inherent sampling 
limitations, it is a reminder that the historical solar occultation-based PSC climatologies are 
strictly representative of their sampling latitudes and not of the polar region as a whole. The 
solar occultation climatologies should be re-evaluated in the context of what we learn from 
CALIPSO. 


Measurement-based estimates of PSC area and volume were derived from the CALIPSO PSC 
observations. The CALIPSO PSC area is about a factor of two smaller than the area with T < 
That, a commonly used proxy to predict PSC occurrence. The CALIPSO PSC area may be 
underestimated since optically-thin PSCs such as NAT haze may not be detected in this initial 
version of the algorithm, but it’s not clear how much these would contribute to the overall 
area. Future versions of the detection algorithm will include polarization measurements to aid 
in identification of this class of PSC. The CALIPSO PSC area is similar in magnitude to the 
area encompassed by T < Tsts. Further study is required to understand the implications of 
this finding and to understand whether T S ts is a better proxy for PSC formation in the 
Antarctic than Tnat or whether the process is too complex for a single parameter proxy to 
capture. In any case, it is clear that CAFIPSO measurements offer a unique opportunity to 
investigate PSC proxies and ultimately develop improved predictors of PSC occurrence. 


The potential of the CAFIPSO data for PSC composition studies was demonstrated. Three 
crude PSC compositions were defined based on combinations of 532-nm scattering ratio and 
volume depolarization ratio measurements. CAFIPSO data can clearly be utilized for 
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discrimination of STS and H 2 O ice PSCs. Perhaps a more important question is the 
composition of the mixed phase clouds. A more in-depth examination of the PSC 
microphysical information imbedded in the CALIPSO data is forthcoming (Poole et al., 
2007). Correlative studies combining CALIPSO and MIPAS PSC observations may provide 
a more comprehensive picture of PSC microphysical properties, particularly the presence of 
large NAT particles in mixed-phase clouds. 
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Table 1. 


CALIOP Instrument Characteristics 


laser: 

wavelengths: 
pulse energy: 
repetition rate: 
receiver telescope: 
polarization: 
footprint/FOV: 
vertical resolution: 
horizontal resolution: 
linear dynamic range: 
data rate: 


Nd: YAG, diode-pumped, Q-switched, frequency 
doubled 

532 nm, 1064 nm 
110 mJoule/channel 
20.25 Hz 
1 .0 m diameter 
532 nm 

100 m / 130 prad 
30-60 m 
333 m 
22 bits 
316 kbps 


Table 2. Spatial Resolution of CALIPSO Lidar Level IB Data Products 


Altitude Range (km) 

Horizontal 
Resolution (km) 

532-nm Vertical 
Resolution (m) 

1064-nm Vertical 
Resolution (m) 

30.1 to 40.0 

5.0 

300 

— 

20.2 to 30.1 

1.67 

180 

180 

8.2 to 20.2 

1.0 

60 

60 

-0.5 to 8.2 

0.33 

30 

60 

-2.0 to -0.5 

0.33 

300 

300 
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Table 3. 


PSC composition definitions. 


Composition 

Scattering Ratio 

Depolarization Ratio (%) 

STS 

<5 

<2 

H 2 0 Ice 

> 10 

> 10 

Mixture 

< 10 

>2 
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Fig 1 . Typical CALIPSO measurement coverage over the Antarctic for a 
single day (28 June 2006). Blue points indicate nighttime orbit segments and 
red points indicate daytime orbit segments. 
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Fraction of Obs. 



Fig 2. Latitudinal distribution of CALIPSO nighttime measurement locations 
for each month of the 2006 Antarctic winter. 
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Fig. 3. Panel (a) shows standard CALIOP backscatter coefficient data from an 
individual orbit on 13 June 2006. Panel (b) shows the same data, but with 
additional smoothing applied to all data below 20.2 km to closely match resolution 
of data above 20.2 km. 
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Temperature, 



Scattering Ratio (532 nm) 

Fig. 4. Scatter diagram of total scattering ratio versus temperature for the 400 - 500 K 
potential temperature layer from an individual orbit of CALIPSO data on 13 June 
2006. 
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Fig. 5. Ensemble distributions of scattering ratio 
June (left panel) and 25 October (right panel). 
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with temperature less than Tpsc- The dashed red line indicates the value of R r . See 
text for discussion. 75-km Averaging 




Fig. 6. Ensemble distributions and scatter diagrams for 25 -km averaging (top 
set of figures) and 75-km averaging (bottom set of figures). See text for 
discussion. 
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Fig. 7. Example of PSC detection for optically thin PSC from 15 June 2006. Top 
panel shows CALIPSO backscatter coefficient image and bottom panel shows the 
corresponding PSC mask produced by the cloud detection algorithm. The level of 
smoothing required to detect the cloud is indicated by the gray scale on the right. 
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Fig. 8. Example of PSC detection for optically thick PSC from 24 July 2006. Top 
panel shows CALIPSO backscatter coefficient image and bottom panel shows the 
corresponding PSC mask produced by the cloud detection algorithm. The level of 
smoothing required to detect the cloud is indicated by the gray scale on the right. 
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Fig. 9. Weekly CALIPSO PSC maps for an altitude of 20 km. The mid-point of each week is 
indicated on individual panels. 
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Fig. 9. Weekly CALIPSO PSC maps for an altitude of 20 km (continued). 
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Fig. 10. Weekly CALIPSO PSC maps for an altitude of 15 km. The mid-point of each week is 
indicated on individual panels. 
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Fig. 10. Weekly CALIPSO PSC maps for an altitude of 15 km (continued). 
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Fig. 11. Longitudinal distribution of PSCs as a function of altitude. Each panel represents weekly 
periods centered on the date shown on the top. See text for discussion. 
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Fig. 11. Longitudinal distribution of PSCs as a function of altitude (continued). 
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Fig. 11. Longitudinal distribution of PSCs as a function of altitude (continued). 
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Fig. 12. Daily time series of CALIPSO PSC frequency as a function of altitude (top 
panel) and corresponding GEOS-4 temperature (bottom panel). 
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Fig. 13. Daily time series of ‘virtual’ solar occultation PSC frequency as a function 
of altitude (top panel) and corresponding GEOS-4 temperature (bottom panel). See 
text for discussion. 
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Fig. 14. Daily time series of CALIPSO-measured PSC area as a function of altitude 
(top panel) and corresponding PSC volume (bottom panel). See text for discussion. 
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Fig. 15. Daily time series of MLS zonal mean nitric acid (top panel), water vapor 
(middle panel), and derived Tnat (bottom panel) for the 5° latitude band centered at 
82.5°S. See text for discussion. 


45 



Altitude, km Altitude km Altitude km 


30 


MLS Nitric Acid Mixing Ratio (ppbv) Mean Latitude = -67 50 



Julian Day 2006 

MLS Water Vapor Mixing Ratio (ppmv) Mean Latitude = -67.50 



160 180 200 220 240 200 280 300 

Julian Day 2006 



Fig. 16. Daily time series of MLS zonal mean nitric acid (top panel), water vapor 
(middle panel), and derived Tnat (bottom panel) for the 5° latitude band centered at 
67.5°S. See text for discussion. 
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Fig. 17. Daily time series of area with T<Tnat using MLS-derived zonal mean profiles 
of Tnat and GEOS-4 temperature data. See text for discussion. 
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Figure 18. Daily time series of area with T<Tsts- See text for discussion. 
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STS Scattering Ratio 


Figure 19. Scatter diagram of CALPSO scattering ratio at 532 nm versus 
volume depolarization ratio at 532 nm for PSC observations from a single orbit 
on 24 July 2006. Regions roughly corresponding to observations of H 2 0 ice 
and STS PSCs are indicated by the blue and red dashed lines, respectively. 
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Figure 20. Examples of PSC composition inferred from the CALIPSO 
measurements. The PSC composition is color coded as indicated by the color 
bar on the right. The top panel corresponds to the tenuous PSC shown in 
Figure 7 and the bottom panel corresponds to the intense PSC shown in Figure 
8 . 


49 


PSC Composition PSC Composition 




